Experiments on the vortex wake of a swimming knifefish 



Zachary J. Taylor^'^, Alexander Liberzon^, Roi Gurka^, Roi Holzman^, Thomas Reesbeck^ 

and F. Javier Diez*^ 

^School of Mechanical Engineering, Tel Aviv University, Tel Aviv, Israel 
^Department of Mechanical and Aerospace Engineering, Rutgers, The State University of 

New Jersey, Piscataway, NJ, USA 
^Department of Mechanical Engineering, Ben-Gurion University, Beer-Sheva, Israel 
^Department of Zoology, Tel Aviv University and the Inter-University Institute for Marine 

Sciences in Eilat, Eilat, Israel 



a: 



C/3 



Abstract 

The knifefish species propels itself by generating a reverse Karman street from an anal fin 
and without significantl y mov i ng its body. This unique feature makes this species' propulsion 
method highly efficient (|BlakeLll983h . It has been suggested that there is an optimal swimming 
range for fish based on the amplitude and frequency of the reverse Karman street. Experiments 
have been performed to measure the ratio between the amplitude and wavelength of vortices 
in the wake of a knifefish. It is suggested that by optimizing the thrust created by the reverse 
Karman street the wave efficiency can be estimated for a given spacing ratio, and present ob- 
servations have an average value of 0.89. The relationship established between spacing ratio 
and wave efficiency, in addition to the measured parameters, will be invaluable for bio-inspired 
designs based on the knifefish. 
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1 Introduction 



In addition to understanding the ecological and evolutionary forces that shape animal diversity, 
the knowledge of fish swimming mechanisms is becoming increasingly important for the design and 
implementation of underwater vehicles. It is natural that designers look to fish for design inspi- 
ration since fish swimming has been ex posed to evolutionary forces that could select for optimal 
swimming (e.g., Sfakiotakis et ah 19991 ) . Fish inhabit a wide array of habitats which place differ- 
ent performance demands. Thus, different fish species evolved swimming behav iours suitable for 
differ ent habitat and task, and there exists a wide range of propulsive methods (jSfakiotakis et all . 
1999l h One of the swim ming ty p es rec ently suggested for use in und erwater vehicles is that k nown 
as gymnotiform motion (Blake, 19831 ). For their swimming robot, Shirgaonkar et al ( 20081 ) have 
taken inspiration from one of the gymnotiform swimmers: the knifefish. In the current study the 
focus is on experiments of a live knifefish species (shown in Fig. [1]). 

The characteristic feature of gymnotiform motion is the anal fin (Fig. [T]) that can be used for 
thrust by creating a travelling sinusoidal wave along the length of its body (Blake, 19831 ). Thus, 
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unlike other forms of fish locomotion, the body can be kept relatively rigid while swimming. While 
this form of swimming does not produce th e greatest th rust, it is known to have high hydrome- 



chanical efficiency ranging upwards of 80% (iBlakel . \l98^ ). In addition, the direction of the wave 
can be reversed implying that the fish can swim backwards as easily as forwards. Knifefish hunt 
for their prey by detecting disturbances in an electric field which suggests one of the reasons for 
this swimming behavior - the prey is detected by swi mming past it and a qui ck change of dir ection 
ensure s that it does not escape after being detected ( Shirgaonkar et al . 20081 ). According to Blake 
(119831 ) it is unclear whether the electric sense preceded the gymnotiform motion in the evolution of 
the species. However, similar swimming behaviours are also observed in other weakly electric fish 
species, suggesting that this swimming mode co-evolved with electrical sense. 

The vortex wake behind fish resembles a Karm an vortex street - only in reverse (jTriantafvllou et all . 

1991 ). Through experiments using a flapping foil, Triantafyllou et al ( 1991 ) suggested that there is a 



narrow range of non- dimensional tail beating frequencie s for optimal thrust between 0.2 < St < 0.4. 
To calculate thrust, IStreitlien and Triantafylloul (119981 ) have reviewed different methods based on 
measurements of the wake and/or body of the fish. One of the three models they used was Karman's 
formulation for the drag induced by a vortex street with the suitable modification of reversing the 
sign of the circulati on. Even though this mode l is re latively simplistic, it was shown to be accu- 
rate to within 10% ( Streitlien and Triantafvllou . 19981 ). One of the challenges of using the Karman 
model experimentally is the apparent need to measure the circulation; however, the thrust of the 
vortex street need not depend on t he circulation directly since the speed of the vortex system. Us, 
and the circulation are linked (e.g., IWeihsl . I1972I ) . If the vortex street thrust coefficient is defined as 
= Ts/^pUja where Tg is the thrust; p the density of water; Uf the speed of the fish; and a the 
vortex waveleng th then the dependency of Ct^ on the spacing ratio b/a can be written similarly to 
BearmanI ( 19671 ) 



Uf) 



^-2) ^coth^ -Fcoth^^ 
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The goal of the current study is to use experimental data to quantify the range of spacing ratios 
that could be expected in the wake of knifefish. These measurements allow the approximation of 
knifefish swimming efficiency, and can h elp establish desig n criteria for those using the knifefish 



species in their bio-inspired designs (e.g.. lCuret et all . I2OIII ). 



2 Details of the experiments 

A sketch of the fish is shown in Fig. [1] identifying the unique anal fin. In order to adhere to 
appropriate policies, regulations and laws for handling and care of vertebrate animals, the tests 
followed protocol 08-017 approved by Rutgers University Animal Care and Facilities Committee for 
studying the African Knifefish (Xenomystus Nigri) swimming methods. Following the protocol, all 
tests were performed only by approved personnel at the Rutgers facility. Tests were performed in 
a recirculating tank constructed out of acrylic that allowed holding the fish in the laser field with 
the possibility of a fiow field to stimulate swimming. A 373 W water pump generated flow speeds 
up to 11 cm/s in the test section with dimensions 8 x 8 x 38 cm^. 

The fiow data were taken underneath the fish in a horizontal plane using two-component Particle 
Image Velocimetry (PIV). The PIV system comprised a Nd:YAG laser (New Wave Research Solo 
PIV 120) and a digital PIV CCD Camera (TSI Powerview Plus 2 MP) at 14.5 Hz. In the current 
study the focus is on six instances when the wake of the knifefish was in the field-of-view (11x8 
cm^) of the camera: three while the fish was swimming forwards (Cases 1-3) and three while the 
fish was swimming backwards (Cases 4-6). It should be noted that although the wake was in the 
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Figure 1: Sketch of a knifefish with the unique anal fin identified. 



field-of-view of the cameras there were only a few snapshots where a sufficient portion of the fish 
appeared in the image to determine its speed. 



3 Results and Discussion 



3.1 Vortex identification 

To measure the spacing ratio in the wake of the knifefish it is necessary to identify the vortex centers. 
Although there is no formal definition of a vortex, locations in the fiow where local pathlines form 
circular patterns are identified by computing the so-called 'swirling strength', Xd. The swirling 
strength is defined as the complex portion of the complex pair of eigenvalues of the two-dimensional 
velocity gradient tensor (if the eigenvalue s are real then ther e is no vortex) , and this method is often 
used on two-dimensional PIV data (e.g., Adrian et al . 200d ) . Once the value of Xd is computed for 
each point in a velocity map, the locations with the highest swirling strength are identified as the 
vortex centers. 

A sample velocity map has been plotted on top of a PIV image with vorticity contours in Fig. 



2(a) and swirling strength contours in Fig. 2(b) Also shown in Fig. 2(b) are the identified locations 
of high swirling strength highlighting the vortex centers. 





Figure 2: An example of an instantaneous PIV image with velocity vectors overlaid in the laboratory 
frame-of-reference. In (a) the contours represent the vorticity field, and in |(b)| are contours of the 
swirling strength, Xd- 



3.2 Vortex Spacing 

One of the classical descriptions of vortex street wakes is the rat i o of w avelength, a, and the cross- 



19671 ). These lengths have been 



stream distance between the two rows of vortices, b (jBearmanl . 
computed according to Fig. [3l The vortex wavelength is computed as the distance between vortices 
1 and 3 and the cross-stream distance, b, is defined as the shortest distance between the vortex at 
(3^2,2/2) and the line connecting vortices at (xi,yi) and {x3,y3). 
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Figure 3: Left: the speed of the knifefish Uf, the speed of the vortices Uw] the width of the body 
d] the numbering of the vortices; the vortex wavelength a; and cross-stream distance between the 
vortex rows h. Right: the analogous circular cylinder moving with the same velocity in the opposite 
direction. 

Table 1: Wake parameters measured for different cases. Each line is the average value for each PIV 
frame. The table is divided into forward (left side) and backward (right side) swimming. 



No. 


b/a 




St 


No. 


h/a 




St 


1 


0.37 


0.79 


0.29 


4 


0.26 


0.90 


0.23 




0.29 


0.87 


0.25 




0.19 


0.97 


0.18 




0.27 


0.89 


0.24 




0.04 


1.00 


0.04 


2 


0.35 


0.81 


0.28 




0.13 


0.99 


0.13 




0.30 


0.86 


0.26 




0.13 


0.99 


0.13 




0.50 


0.72 


0.36 




0.08 


1.00 


0.08 




0.22 


0.94 


0.21 


5 


0.17 


0.98 


0.17 




0.44 


0.74 


0.33 




0.32 


0.84 


0.27 




0.25 


0.91 


0.23 


6 


0.14 


0.99 


0.14 




0.19 


0.96 


0.18 












0.26 


0.91 


0.23 










3 


0.47 


0.73 


0.34 












0.26 


0.91 


0.23 












0.41 


0.76 


0.31 












0.16 


0.98 


0.16 












0.24 


0.92 


0.22 











The number of vortex triplets used to estimate the spacing ratio ranged between 1-4 in each 
PIV frame, and each frame-average spacing ratio is presented in Table [TJ A significant range of 
spacing ratios i s observed in Ta ble [1] which agrees with the variety of spacing ratios reported in the 
literature (e.g., Bearman . 19671 ). One of the particular features of the knifefish is that the vortices 



are created by the sinusoidal shape and movement of the fin rather than the beating of its tail. Thus, 
rather than exclusively creating vortices behind its body, the knifefish creates vortices underneath 
its body (e.g.. Fig. [2]). This swimming means it can instantaneously change the spaci ng ratio of 
the vo rtices since it is governed by the shape of the anal fin. For fish and robotic fish, lEpps et al 
(j2009l ) reported a wide range of spacing ratios (0.2-0.8) which range noticeably higher than those 
found for the knifefish. 



3.3 Efficiency 

The research on vortex street wakes by Bearman ( 19671 ) used Karman's drag equation for a vortex 



street and the BCds/ d{h/ a) = stability criterion proposed by Kronauer to link the spacing ratio 
and the convection speed of vortices in the wake of bluff bodies. The schematic shown in Fig. 
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[3] shows the reverse Karman street generated in the wake of the knifefish and the corresponding 
motion in the opposite direction of a circular cyhnder. Thus, if the fish generates the same vortex 
arrangement as the cyhnder while moving in the opposite direction, the drag induced by the vortex 
street is analogous to the thrust produced by the fish (Eq. [1]). Moreover, re- framing the Kronauer 
stability criterion to maximize thrust (instead of minimizing drag) based on the vortex sp acing then 
the sp eed of the vortex system, Ug, can be related to the speed of the fish, Uf, similarly to iBearman 
(jl967l ). Differentiating Eq. [T]with respect to the spacing ratio b/a and setting it to zero yields Eq. 
[2] where k = nb/ a is substituted for compactness. 



Uf 

u. 



Tf = 2 1 + 



sinh^A; — ktanhk 



One of the classical measures of swimming efficiency is the so-called 'wave slip', or wave efficiency, 
defined as the ratio Uf /U,„ where is the wave speed, or convection speed, of the vortex street 
( Sfakiotakis et al , 19991 ). In the assumptions behind Eq. [H the wave speed of the vortices is defined 
as Uw = Uf + Us- Thus, the wave efficiency is calculated using Eq. [2] and 



El 



1 + Us/Uf 



(3) 



These data are added to Table [T] for each PIV measurement in al l 6 cases. Th e wave efficiency is not 
the same as the hydromechanical efficiency cited previously by Blake ( 19831 ) . However, t he ra nge 



of wave efficiencies found here are roughly consistent with those reported by ICuret et all (j201ll ) of 
0.6-0.9. The efficiencies approach unity as the spacing ratio approaches zero. However, it is not 
possible for a fish to continue swimming if 6/a — ?> 0. As previously mentioned, since the fish did 
not often appear in the image at the same time as its wake, it was not possible to estimate Uf 
very often from the current dataset. However, when it was possible, the wave slip was computed 
using the convection speed of the vortices as the wave speed (taken as the speed at the center of the 
vortex). The measured wave efficiency ranged between 0.85-0.89 which agrees with the predicted 
values based on the spacing ratios in Table [H 

In addition to the wave slip , it is well known that th ere is an optimal Strouhal number (St = 
fb/Uf) range of 0.2 < St < 0.4 (jTriantafyllou et all . ll99ll ) for propulsive wakes. Since the frequency 
of the vortex street is / = U^/a the Strouhal number can be written as 



St 



Uu^b 
Ufa 



(4) 



The Strouhal numbers have been added to Table [T] using the measured spacing ratios and Eq. 
51 and for most cases th e calculated Strouhal numbers fall inside the optimal range predicted by 



Triantafyllou et "al ( 1991 ). The cases falling outside of this optimal range correspond to cases when 
b/a 0. 



4 Conclusions 

Through PIV experiments on a swimming knifefish the spacing ratio of the wake is calculated for 
several repeated cases. The PIV measurements also show that the knifefish can create a thrust- 
generating vortex street wake while a portion of that wake remains underneath its body. This ability 
is due to the unique anal fin and one of the reasons for its superior swimming efficiency. A link 
has been suggested between spacing ratio and efficiency (Eq. [3|), thus it is not surprisi ng that the 
spacing ratios measured for the knifefish are slightly smaller than in previous studies (jEpps et all . 
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20091 ) since knifefish are known to be amone; the most efficient swimmers (jBlakel . ll983l ). While the 



two-dimensional ap proximations seem to rn atch previously published wave efficiencies and optimal 
Strouhal numbers, Shirgaonkar et al ( 20081 ) have shown the importance of the three-dimensional 
connections between the vortices and further studies are required to address how the neglected 
three-dimensionality affects the swimming efficiency of knifefish. 
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